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Abstract 

By extending the statistical distributions to the transverse degree of freed- 
hom, we account for a multiplicative factor in the Fermi-Dirac functions of 
the light quarks, we were led to introduce in a previous work to comply with 
experiment. We can also get light antiquark distributions, similar to those 
we proposed earlier. 
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Resulting from a research which begun many years ago PP, we proposed 
for the parton distributions some expressions inspired by quantum statis- 
tics, with a relationship between the light q and q distributions following from 
the chiral properties of QCD |3j| . We determined the few free parameters in- 
volved, from a selected set of precise deep inelastic scattering data [21 • Our 
approach has a strong predictive power, because, once the parameters, occur- 
ing in the better known unpolarized u{x) and d{x) distributions, are fixed, 
the polarized distributions, /S.u{x) and A(i(a;), as well as the antiquark dis- 
tributions u{x), d{x), Au{x) and Ad{x), are also fixed. 
First, we made a prediction for a positive Ai{x, Q"^) at high x jl], which was 
successfully verified by recent data at Jefferson Lab [3]. Next, we predicted 
for the high x behavior of the ratio d{x)/u{x) a monotonously increasing 
function of x. An earlier result from E866 at FNAL [B], seemed to contradict 
our prediction, but the comparison with the most recent results on Drell-Yan 
processes by the NUSEA [7j, shows that this is not the case jH]. Finally, the 
agreement with the data on e^p neutral and charge current reactions [H] PU] ) 
is particularly successful. In Ref. [HI we stressed the experimental evidence 
in several structure functions for a typical property of the statistical distri- 
butions: the change of slope for x > Xq^, where X^^ denotes the largest 
"thermodynamical potential" , whose definition and value are given below. 
We now recall some of the basic features of the statistical approach. The 
fermion distributions are given by the sum of two terms P|, a quasi Fermi- 
Dirac function and a helicity independent diffractive contribution equal for 
all light quarks: 

^2. ^^oV Ax' 

exp[(x — AQg)/xJ + 1 exp(a;/x) + 1 



V n2^ _ MXp^T^x'' Ax' 
xq {X, Qo) - ^^p^^^ ^ ^ ^ + ^^p(^/-) ^ ^ , (2) 



at the input energy scale Ql = 4GeV^. 

The parameter x plays the role of a universal temperature and X^^ are the 
two thermodynamical potentials of the quark q, with helicity h = ±. For the 
gluons we consider the black-body inspired expression 



xG{x,Ql)= ff'' , (3) 
exp(x/x) — 1 



2 



a quasi Bose- Einstein function with be = b + 1, since we believe that the 
diffractive contribution in Eqs. (1,2) is strongly related to the gluons. 
We also assume that, at the input energy scale, the polarized gluon distribu- 
tion vanishes, so 

xAG{x,Ql) = 0. (4) 

For strange quarks we assumed that s = s, but we did not introduce any 
more parameters and we simply related them to u + d to comply with ex- 
perimental data PI]. Similarly we took As = As and we related them to 
{Ad — Am), to agree with the second Bjorken sum rule. 
The eight free parameters in Eqs. (1,2) {A, A and Aq are fixed by the nor- 
malization conditions u — u = 2, d — d = l and by the proton momentum 
sum rule) were determined at the input scale from the comparison with a 
selected set of very precise unpolarized and polarized DIS data from NMC, 
BCDMS, E665, ZEUS, HI for F^''^{x,Q^), from CCFR for xF^^{x,Q^) and 
from SMC, E154 and E155 for ^f'^'"(x,Q2). 
We obtained 

X = 0.09907, b = 0.40962, b = -0.25347 , (5) 



= 0.46128, Xq- = 0.29766, X^" = 0.30174, = 0.22775 , (6) 

A = 0.08318, A = 1.74938, A = 1.90801, Aq = 14.27535 . (7) 

In Eqs. (1,2) the factors Xq^ and {Xq^^)~^ introduced in the numerators of 
the g's and g's distributions, imply a modification of the quantum statistical 
form, we were led to introduce to agree with experimental data. As we will 
see now, the justification of these factors can be understood by studying the 
role of the transverse degree of freedhom for the partons. 
Let us denote Pi{x,p'^) ^ a parton (quark or gluon) distribution with trans- 
verse momentum pj-, which is assumed to be much smaller than its longitu- 
dinal momentum xpz, where Pz is the momentum of the proton with mass 
M. We first consider the momentum sum rule, which reads now 

^ dx xpi{x,p^)dp% = 1 , (8) 

j ^0 ^0 

^Here i stands for the quark flavour and for the spin component of the partons; pi (a;, p^^) 
has now the dimension of GeV~^. 
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where the upper hmit for the pip integral, {j)'ip)maxi will be determined below. 
Next, we consider, the condition for the energy sum rule, analogous to Eq. (8), 
where we assume that the proton mass is much smaller than its momentum, 
an approximation which holds in the deep inelastic regime. Consequently, 
the proton energy is Pz + M"^ /2pz., similarly the energy of a massless parton 
is xpz + p^/2xpz, so this energy sum rule reads 

/•I j'{p\<)max 

dx j xpi{x,PT)dpT + 

j Jo Jo 

\ Jo Jo ^^^^'^^^^^^^ = + ^ • 

By using Eq. (8) in Eq. (9), it is easy to derive the constraint 

2_, dx pi{x,pl)^dpl = M'^ . (10) 

^ Jo Jo ^ 

By comparing Eq. (10) and Eq. (8), where it is clear that x is bounded by 1, 
we see that p^/x is bounded by M^, so {p'^)max = xM"^ . 
The above sum rule will be used to determine the pt dependence, as we will 
explain now. We recall the general method of statistical thermodynamics 
[12j . to find the most probable occupation numbers nj, for the energy levels 
ej, when the total energy of distinguishable particles is E. Since we have 

N = '^^ni and i?=^^njej, (11) 

i i 

one should look for the maximum of 

ln[iV!/ Y[n^]+a{N-Y,n^)+ (3{E - n^e^) , (12) 

ill 

with the Lagrange multipliers method and we find 

Hi = exp(-a - f3ei) . (13) 

The Lagrange multipliers a and j3 are fixed by the two above constraints 
Eq. (11). Now by putting in correspondence rii with pi{x,p'^), N with 1 and 
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a with for the sum rule (8), and E with M^, with p^/x and /? with 
for the sum rule (10), one finds, in correspondence with Eq. (13) 

2 

Vi{x,vV) = exp(— + ^) . (14) 

is a parameter which has the dimension of a mass and it will be determined 
later. We recall that Eq. (14) describes well the pt behavior of the particles 
produced in very high energy hadron-hadron scattering [13j. For fermions, 
the Boltzmann exponential form Eq. (14) is expected to be modified into the 
product of two Fermi-Dirac expressions. 

For the sake of clarity, we will first ignore the pt dependence of the diffractive 
contribution of quarks, antiquarks (second terms in Eqs. (1,2)) and of the 
gluons, which dominate the very low x region. It certainly corresponds to 
very small pr values, since {pT)max = xM^ as seen above. So let us consider 
now the non-diffractive terms and instead of the helicity, as in Eqs. (1,2), the 
spin component of the partons along the momentum of the proton, which co- 
incides with helicity only when pt = 0. The quantum statistics distributions 
for quarks and antiquarks read in this case 

S / 2 \ (^CC^ 1 , , 

xq ^{x,p^) = ^^^^^^^ _ ^ [exp[(p2 - Y,f)/x] + 1] ' ^^^^ 



[exp[(x + Xo-/^)/x] + 1] [exp[(p2 /a;/i2 + Y,-/^)/x] + 1] ' 



xf^{x,pl)= ^ , ^ , , (16) 



where F{x) and F{x) do not depend on spin and flavour, X and Y are the 
thermodynamical potentials related to the sum rules (8) and (10), respec- 
tively ^. At high Pt, one has a Gaussian behavior, with a width proportional 
to ^/x. It is important to note that this is at variance with the usual fator- 
ization assumption of the dependences in x and pt [IH CS] • We observe that 
the upper limit of the integral in x, means that one parton is taking all the 
proton momentum. Similarly {pT)max = xM"^ for the p\ integral corresponds 
to the situation where one parton is taking all the transverse energy. There- 
fore to simplify the formulas, we shall take as upper limit infinity, since the 

''Here for convenience, we have divided by x the argument of the exponential in the 
Fermi-Dirac expression of the pT dependence, in agreement with the expression for the x 
dependence. 
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contribution of the tail is negligible. In fact when one now integrates over 
p^, the right-hand side of Eq. (15), one finds 



exp[(p2./a;/i2 - Y^^„^)/x)] + 1 



, , - -x/x^xLii(-exp[ro^;/x]) . (17) 



Here Lii denotes the polylogarithm function of order 1, which is known to 
arise from the integral of Fermi-Dirac distributions and is such that 

-Lii(-e^)= r^^J^^=ln{l + ey). (18) 

It is reasonable to assume the proportionality relationship 

Yof = kXl; . (19) 

It implies that the partons with a larger contribution to their first moments 
from the non-diffractive part, not only have a broader x dependence, but also 
have a broader dependence, at every x. We recall that Lii is an increasing 
function of its argument and we see from Eq. (18) that for large values of 
?/, it becomes approximately proportional to y and one has the high degen- 
eracy of the Fermi gas. Therefore this fully justifies, the phenomenological 
assumption of the proportionality to Xq^ made in Ref. j2], as seen in Eq. (1). 
By taking in Eq. (15) ^ 

Ui{-exp[Y^l/x\)fi^x 

we recover the first term in Eq. (1) for the quark, which is the dominating 
parton at large x. Eq. (20) implies for the other quarks, that instead of Xq^, 
one should use the factor X^Lii(— exp[FQg/a;])/Lii(— exp[Y'Q;J^/a;]). 
As long as the non-diffractive term for the antiquarks, by integrating on p'^, 
one finds the factor Lii(— exp[— Fq^/x]). At large negative y, Lii becomes 
proportional to exp(— y), which corresponds to the Boltzman limit. Like 
for the quarks, we will fixe the normalization using the largest distribution, 
namely d~, so similarly to Eq. (20) we will take in Eq. (16) 

= • (21) 

Uii-expl-Y^l/xDfi^x 



*We identify X^^ with Xl}^ given in Eq. (6) and similarly for the potentials Y^^' . 
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We recover the first term in Eq. (2) for d and for the other antiquarks instead 
of {X^^)-\ one should use (X+)-iLii(- exp[-yoV^])/Lii(- exp[-r+/x]). 
Now, in order to determine /c, we require no change for the quark distribu- 
tions, once integrated over p^. This can be achieved by choosing 

k = 1.42 . (22) 

This value indicates that the potentials Yq^ are of the same order as Xq^ and 
moreover it implies lower values for the non-diffractive contributions of , 
and u" . Since the antiquark distributions are dominated by the diffractive 
contribution, these small modifications are still consistent with the data. 
At this stage we can try to use the energy sum rule Eq. (10) to determine /i, 
the only unknown parameter, so far. We find to a very good approximation, 
the remarkable value, 

H = IGeV . (23) 

Moreover the energy sum rule is found to be saturated by a quark contribu- 
tion of 98%, whereas the antiquark contribution is only 2%. This value of ix 
was found by neglecting the gluon and the diffractive contributions and it is 
clear that if they turn out to be important, fi will be smaller. 
For completeness we now turn to the universal diffractive contribution to 
quarks and antiquarks in Eqs. (1,2), namely xq^{x, Qq) — Ax^/[exp{x/x)+l\. 
Since b < (see Eq. (5)), the introduction of the pt dependence cannot be 
done similarly to the non diffractive contributions, because in the energy sum 
rule Eq. (10), it generates a singular behavior when a; — > 0. Therefore we 
propose to modify our prescription by taking at the input energy scale 

D 2 Ax''^'^ 1 1 

{x,Pt) = in(2);,2^ [exp(x/x) + 1] [exp(p2 /a;a;2^2) + ^ ' (24) 

whose Pt fall off is stronger, because x//^ is now replaced by x'^/j,'^. Note that 
this is properly normalized to recover a;g^(a;, Qq) after integration over py. 
We have checked that xq^{x^p'^) gives a negligible contribution to Eq. (10), 
as expected, and we are led to a similar conclusion for the gluon, which 
has the same small x behavior. However since the gluon is parametrized 
by a quasi Bose-Einstein function, one has to introduce a non-zero potential 
Yg, to avoid the singular behavior of Lii(exp[— 1^/2^]), when Yq — Q (see 
Eq. (18)). The value of Yq is not constrained, but by taking a very small 
Yg, it does not affect the energy sum rule. Of course, one may ask why the 
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gluon contributes so much to the momentum sum rule, Eq. (8), and so httle 
to the energy sum rule, Eq. (10). Using a better regularization procedure, 
we might find a larger contribution and consequently a smaller value. For 
this reason, we give the pt dependence as a function of the dimensionless 
variable p^/fi"^. 

For illustration, we display in Fig. 1, xjj?u{x^p'^/ jj? ^Ql), the predicted di- 
mensionless u-quark parton distribution at input energy scale, versus p%/ fx'^ 
for different x values. Note that the curves are limited to the value {p'^)max = 
xM"^ and the Gaussian behavior in the high px region follows a rather flat 
behavior for 10~^ < < 10~^, or so. This feature is a natural conse- 

quence of the statistical approach and it coincides, for a mean x value, with 
the simplifying assumption made in recent papers (see e.g. [13]) without real 
justiflcation. For p%/ < 10~^, the very low x region is strongly enhanced by 
the diffractive contribution, due to the factor x^~'^ (see Eq. (24)). However, 
since the regularization procedure proposed above is not unique, it should 
not be considered as a deflnite prediction. 

The theoretical fundation of the statistical parton distributions proposed 
earlier in Ref. [2] is strengthened by the present extension to the transverse 
degrees of freedhom. Indeed the energy sum rule Eq. (10) cuts naturally 
high p]. values and this is very appealing, because an increasing phase-space 
for p^ linear in Q"^ would support the absence of quantum statistical effects 
for diluted fermions and bosons. Our predicted light partons distributions 
Pi^x^pip), both unpolarized and polarized given by Eqs. (15,16), remain to 
be checked against appropriate experimental data, sensitive to some effects 
arising from the transverse momentum of partons inside the proton. Here it 
is important to emphasize that the pt dependence is only valid in a limited 
kinematic region, say 0.1 < a; < 0.9, which is not dominated by the diffractive 
contribution. Clearly in this large x region, we predict the mean value of 

namely < >, to increase with x, a behaviour which is at variance 
with what one expects in the very small x region from the BFKL evolution, 
assuming the emission of a hudge number of gluons JHl- However, so far 
there is no experimental evidence for the BFKL evolution. 

Finally, we note that our procedure restricts to rotational invariance, the 
transverse momentum dependence of the parton distributions. The next step 
should be to introduce some azimuthal dependence, which is needed to study 
the important topic of single-spin asymmetries |17j . 
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Figure 1: The predicted xfi'^u{x,p'^/ ii'^,Ql) statistical distribution, versus 
Px/l^''^ for different values of x. 
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